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ABSTRACT 

Aims. The Strouhal number (St), which is a nondimensional measure of the correlation time, is determined from numerical 
models of convection. The Strouhal number arises in the mean-field theories of angular momentum transport and magnetic field 
generation, where its value determines the validity of certain widely used approximations, such as the first order smoothing 
(hereafter FOSA). More specifically, the relevant transport coefficients can be calculated by means of a cumulative series 
expansion if St < St cr it ~ 1. 

Methods. We define the Strouhal number as the ratio of the correlation and turnover times, which we determine separately, the 
former from the autocorrelation of velocity, and the latter by following test particles embedded in the flow. 
Results. We find that the Strouhal numbers are, generally, of the order of 0.1 to 0.4 which is close to the critical value above 
which deviations from FOSA become significant. Increasing the rotational influence tends to shorten both timescales in such a 
manner that St decreases. However, we do not find a clear trend as function of the Rayleigh number for the parameter range 
explored in the present study. 



Key words, convection - hydrodynamics - dynamo theory 



1. Introduction 

The mean-field theories of angular momentum transport 
(e.g. Riidiger 119801 I1989f) and hydromagnetic dynamos 
(e.g. Steenbeck & Krause ITM?! Krause & Radler IrMTjl 
require knowledge of the Reynolds stresses and the mean 
electromotive force, respectively. Direct analytical calcu- 
lation of these quantities for astrophysical purposes is not 
possible at present due to the lack of an established the- 
ory of turbulence. Numerical calculations are not much 
better off since in most astrophysical systems the compu- 
tational resources needed in order to resolve all physically 
relevant scales are several orders of magnitude larger than 
currently available. 

In mean-field theories these problems are circumvented 
by relating the Reynolds stresses and electromotive force 
to the mean quantities (the rotation vector 17 and the 
mean magnetic field (B), respectively) and their gradi- 
ents by means of a cumulative series expansion (see van 



Kampen [T9"74afbT) . In the best known and most often 
used approach, the first order smoothing approximation 
(FOSA), only the first terms of these expansions are taken 
into account. This approach can be shown to be valid if 
either the relevant Reynolds number, Re = ul/v , Rm = 
ul/rj, or the Strouhal number 



= \(U • V)u\ _ ^ ^ Tc_ 

\du/dt\ ~ l c t to 



(1) 
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is sufficiently small. Above, u and / are the typical veloc- 
ity and length scales, and r c and l c are the correlation 
time and length of the turbulence, respectively. We show 
in Sect. 13. 31 that l/t to oc u/l c . Only the case of small St 
is relevant in astrophysical circumstances, where typically 
Re, Rm ^> 1. Hence the condition St <C 1 should be ful- 
filled if FOSA is to be valid in stellar convection zones. If 
this is not the case then FOSA is likely to be too crude a 
truncation but, as long as St is smaller than the critical 
value for convergence it is possible to construct higher- 
order mean-field theories by including more terms of the 
cumulant expansion. However, a universal critical value 
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of St for convergence of the cumulant expansion does not 
exist because this depends on the geometry of the prob- 
lem and on the flow itself. In general though, the critical 
value is thought to be less than unity. For simple turbu- 
lence models it is possible to compute it; Nicklaus & Stix 
(1988 , hereafter NS88) obtained Stent = 1 for their model. 
Furthermore, the results of NS88 show that the higher or- 
der effects remain small in comparison to the FOSA result 
when St < 0.5. 

It is well-known that the requirement for a small 
Reynolds number is not satisfied in stellar environments. 
However, the question of the Strouhal number is not set- 
tled. On account of the observations of the solar surface 
granulation, one can estimate the correlation and turnover 
times to be roughly equal, indicating that St « 1 at the 
solar surface (see e.g. Chapter 6 of Stix 2002). Similar val- 
ues can also be estimated for supergranulation for which 
typical numbers are u ~ 100 m s" 1 , I w 10 7 m, and 
t 10 5 s. However, even for granulation the precise value 
has, to our knowledge, not been established, and noth- 
ing is known about St in the deeper layers. Furthermore, 
recent results from forced turbulence calculations indi- 
cate that if the higher order correlations in the equa- 
tions of the passive scalar flux (Brandenburg et al.|J004 
hereafter BKM) and electromotive force (Brandenburg & 
Subramanian I2()()5|l are taken into account via the so- 
called minimal r-approximation (Blackman & Field 2002 
2003), the Strouhal number can be seen to substantially 
exceed unity and roughly equal to unity, respectively. 

Although forced turbulence is rather different in com- 
parison to convection, the aforementioned studies still 
raise the question whether the results of the convection 
calculations can be interpreted within the framework of 
the standard mean-field theory as has been done in nu- 
merous studies during recent years (e.g. Brandenburg et 
al. ITWil Pulkkinen et al. IHMl Ossendrijver et al. I27M1 
l2T)(T2l Kapyla et al. l2TMl Rudiger et al. l2T)0"5)> . Motivated 
by the unknown status of the Strouhal number for convec- 
tion and the previous studies on the subject in different 
contexts, we set out to calculate St from numerical calcu- 
lations of convection. In order to do this, we calculate the 
correlation time from the autocorrelation of velocity and 
determine the turnover time by following test particles 
embedded into the flow and define the Strouhal number 
as the ratio of the two. 

The remainder of the paper is organised as follows: 
Sect. |5] summarises briefly the numerical model used and 
in Sect.|31the results of the study are discussed. Finally, 
Sect. 0] gives the conclusions. 

2. The model 

A detailed description of the convection model can be 
found in Kapyla et al. I|2004l hereafter Paper I) and the 
convection calculations are made with a setup identical to 
that used in Paper I. See Table ^ for a summary of the 
main parameters. 



Table 1. Summary of the calculations and the main 
parameters. From left to right: the Rayleigh, Reynolds, 
Taylor, and Coriolis numbers, the latitude, and the grid 
size. 



Run 


Ra 


Re 


Ta 


Co 




Grid 


ICoO 


1.25 


10 b 


95 








— 


T77\ 

48 3 


CoO 


2.5 ■ 


10 5 


140 










64 3 


mCoO 


5.0 ■ 


10 s 


190 










96 3 


hCoO 


10 6 


246 










128 3 


CoOl-00 


2.5 ■ 


10 b 


140 


203 


0.10 


0° 


64 3 


CoOl-30 


2.5 • 


10 5 


140 


203 


0.10 


-30° 


64 3 


CoOl-60 


2.5 ■ 


10 s 


138 


203 


0.10 


-60° 


64 3 


CoOl-90 


2.5 ■ 


10 5 


138 


203 


0.10 


-90° 


64 3 


Col-00 


2.5 ■ 


10 5 


139 


2.03 ■ 10 4 


1.04 


0° 


64 3 


Col-30 


2.5 ■ 


10 5 


145 


2.03 • 10 4 


1.00 


-30° 


64 3 


Col-60 


2.5 ■ 


10 5 


141 


2.03 • 10 4 


1.03 


-60° 


64 3 


Col-90 


2.5 ■ 


10 5 


139 


2.03 • 10 4 


1.05 


-90° 


64 3 


ColO-00 


2.5 ■ 


10 5 


337 


2.03 • 10 6 


4.24 


0° 


96 2 x 64 


ColO-30 


2.5 ■ 


10 5 


121 


2.03 ■ 10 6 


11.8 


-30° 


96 2 x 64 


ColO-60 


2.5 ■ 


10 5 


105 


2.03 • 10 6 


13.6 


-60° 


96 2 x 64 


ColO-90 


2.5 • 


10 5 


104 


2.03 • 10 6 


13.7 


-90° 


96 2 x 64 



For the purposes of the present study we have added 
the possibility to follow the trajectories of Lagrangian test 
particles in the model. In order to integrate the trajectory 
we need the velocity at the position of the test particle 
at each integration step. This is done by finding the grid 
points (n x , n y ,n z ) next to the test particle and using linear 
interpolation to obtain the velocity at the correct position 



u(x tp ) = u(n) - ^2 7T- U ( n ~ e *) ' ( 2 ) 

i—x.y,z 

where x tp — (xt p , ytp, z tp ) is the position vector of the 
test particle, n = (n x ,n y ,n z ) denotes the grid point 
next to the test particle, e.{ the unit vector in direction 
i, 5xi = Xi — Xtp the distance between the test particle 
and the grid point next to it, and Axi the grid spacing 
in direction i. In the present calculations we follow one 
thousand test particles that are introduced at the middle 
of the convectively unstable layer at random horizontal 
positions after convection has reached a statistically sta- 
tionary state. 

The calculations were made with a modified version 
of the numerical method described in Caunt & Korpi 
(2001]). The calculations were carried out on the KABUL 
and BAGDAD Beowulf clusters at the Kiepenheuer- 
Institut fur Sonnenphysik, Freiburg, Germany, and on the 
IBM eServer Cluster 1600 supercomputer hosted by CSC 
Scientific Computing Ltd., in Espoo, Finland. 
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Fig. 1. Horizontally averaged correlations of u z with re- 
spect to the snapshot at t = 100 for the run CoO. Each 
curve is separated by At = 1 in time units of \Jdjg. 

3. Results 

3.1. Correlation time 

We determine the correlation time from the velocity au- 
tocorrelation function 



C[ui{x,t ),Ui{x,i)} 



(Uj(x,t )Ui(x,t)} 

^(uj(x,t Q ))(uj(x,t)) ' 



(3) 



where i denotes the velocity component used, to and t 
are the times from which the snapshots were taken, and 
the brackets denote horizontal averaging. We estimate the 
correlation time, r c , to be the time after which the corre- 
lation drops below a fixed threshold value, in this case 0.5. 
However, the correlation time still depends on depth, and 
there are discrete time intervals between the stored snap- 
shots (see Figure P). To remedy the latter, we calculate 
the correlation time within the convection zone for each 
depth and use linear interpolation to find a more accurate 
value for the time when C — 0.5. Furthermore, it makes 
sense to average the correlation time over the convectively 
unstable layer since convection is more of a global rather 
than local nature with the present parameters (see Figure 
|3J). To check the time dependence of r c we calculate it 
with respect to about two hundred snapshots for each cal- 
culation. The final correlation time is an average of these 
individual values. 

Figure ^gi ves an example from run CoO, showing the 
horizontally averaged correlations of the vertical velocity 
u z from twelve snapshots, each separated by one time unit, 
with respect to the snapshot at to = 100. The correlation 
diminishes monotonically as a function of time and for the 
eleventh snapshot the correlation is below 0.5 in the whole 
convection zone. Using the procedure described above, we 
find the correlation time to be r c « 9.8 for this snapshot 
(average over the correlation times with respect to 206 
different snapshots gives a value 9.4, see the second col- 
umn of Table 01 ■ If the correlation is calculated for one of 



the horizontal velocity components, t c is similar in magni- 
tude, except near the surface and in a layer immediately 
below the convection zone where it is somewhat longer 
than the one calculated from the vertical velocity. This 
effect can be understood to arise from the persistent hor- 
izontal flows near the boundaries of the convection zone 
where the up- and downflows diverge to the horizontal di- 
rections. In what follows we estimate the correlation time 
from the vertical velocity, for which r c remains more or 
less constant within the convection zone as indicated by 
Figure [I] As stated above, our final result is an average 
over snapshots and the corresponding standard deviations 
of the correlation times are given in Table 

3.2. On the correlation length 

In principle, the correlation length of the turbulence can 
be determined from Eq. J2J using r = |xi — x%\ as the ar- 
gument instead of the time. In order to do this, we choose 
one thousand random grid points within the convectively 
unstable layer and compute the crosscorrelations of the 
velocities between all the points. We do this procedure for 
the same snapshots from which the correlation times were 
calculated, and average over all snapshots. Due to the fi- 
nite spatial resolution of the calculations the smallest dis- 
tance between two points cannot be smaller than the grid 
spacing. Thus the correlations are binned, each bin cover- 
ing a range Ar = 0.075d. A characteristic result is shown 
in Figure|21 where the correlations from runs with different 
Coriolis numbers at the southern pole are presented. The 
main conclusion is that the correlation diminishes faster 
as function of r when rotation is more rapid, implying that 
l c decreases as a function of rotation. There is, however, 
a problem of how exactly to define l c . Due to the bin- 
ning the average correlation in the bin with the smallest r 
varies from one calculation to the other, and setting a fixed 
threshold value would be more arbitrary than in the case 
of the correlation time. We see that a reasonable solution 
to bypass this problem is to use test particles in order to 
determine the turnover time, tt (ct l c /ut), that captures 
the changing spatial scale of convection as a function of 
rotation. 

3.3. Turnover time 

The turnover time is usually estimated by dividing a char- 
acteristic length scale by a characteristic velocity. For con- 
vection calculations one often uses the depth of the con- 
vection zone divided by the volume-averaged rms velocity, 
i.e. 



M _ 



d_ 



(4) 



where the superscript s refers to 'simple estimate'. For ex- 
ample, for the CoO run, u% « 0.085 and d = 1, giving 

(s) 

tlJ ss 12.0. Using this estimate with r c = 9.4 would in- 
dicate that St is approximately 0.8. Although this simple 
method probably gives the correct order of magnitude of 
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Fig. 3. Vertical positions of two test particles in run the 
CoO. The dotted vertical lines denote the times where the 
particle denoted by the solid line changes direction from 
downward motion to upward motion. The dash-dotted line 
shows a particle which is stuck in the lower overshoot layer 
for the whole duration of the calculation. 



the turnover time, it is still quite a crude estimate since 
the values for the typical scales and velocities are rather 
uncertain and vary nonlinearly as function of rotation and 
Rayleigh number (see below) . Thus a more precise way of 
calculation is desirable. 

A way to improve the estimate given by Eq. is to 
follow the trajectories of Lagrangian test particles in the 
flow. This is done by finding the times where a particle 
changes its direction, i.e. turns over. Thus, one turnover 
time would be the time between two consecutive turns 
to the same direction, e.g. from downward motion to up- 
ward motion. The advantage of this method is that the 
assumption of the vertical scale of convection is removed. 
However, the danger with this method is that the small- 
est scales begin to dominate due to, for example, contribu- 



Table 2. From left to right: rms-velocity averaged over 
the convectively unstable region and time, correlation time 
from the autocorrelation of the vertical velocity, turnover 
time from test particle trajectories, and the Strouhal num- 
ber. The last column states the number of snapshots with 
respect to which the correlation times and lengths were 
calculated. 



Run 






t to 


St 


N 


ICoO 


0.085 


11.2 ± 2.3 


25.9 


0.43 


209 


CoO 


0.084 


9.4 ±2.1 


25.1 


0.37 


206 


mCoO 


0.080 


8.8 ±2.6 


24.0 


0.37 


212 


hCoO 


0.073 


8.9 ±1.6 


23.4 


0.38 


199 


CoOl-00 


0.084 


9.4 ±2.2 


25.0 


0.38 


206 


CoOl-30 


0.084 


9.2 ± 1.5 


25.1 


0.37 


202 


CoOl-60 


0.082 


9.8 ± 2.9 


25.0 


0.39 


207 


CoOl-90 


0.082 


9.0 ± 1.6 


25.5 


0.35 


203 


Col-00 


0.083 


7.3 ± 1.4 


24.5 


0.30 


207 


Col-30 


0.086 


6.1 ±0.9 


23.4 


0.26 


198 


Col-60 


0.084 


6.9 ±0.8 


23.7 


0.29 


207 


Col-90 


0.083 


7.1 ±0.9 


23.8 


0.30 


206 


ColO-00 


0.201 


2.2 ± 0.3 


12.6 


0.17 


202 


ColO-30 


0.072 


3.9 ±0.3 


21.0 


0.19 


200 


ColO-60 


0.062 


3.4 ±0.2 


20.5 


0.17 


201 


ColO-90 


0.062 


3.1 ± 0.4 


18.4 


0.17 


201 



tions from particles stuck in the stable layer. This problem 
is discussed below. 

Figure |31 shows the trajectories of two test particles 
in the run CoO. The solid line shows a particle which is 
carried by the convective flow throughout the calculation 
(the particles were introduced in the flow at t = 60). In 
the figure we only denote the turnovers from downward to 
upward motion, but including also the opposite changes 
of direction a total of 21 turnovers are registered. Figure 
0] shows a histogram of the registered turnover times from 
run CoO. In total, the thousand particles make 21300 in- 
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dividual turnovers in the course of the calculation. The 
distribution is centered near 20 having an exponential tail 
towards longer times. Arithmetic average over the distri- 
butions shown in Figure 0] gives a turnover time ito ~ 25, 
which indicates that the Strouhal number to be signifi- 
cantly less than 0.8 which was obtained with the simple 
estimate of the turnover time. 

The test particle method takes into account the vari- 
able spatial scale of convection but it also picks up the 
small scale turnovers in the stable layer. A small number 
of particles get stuck in the lower overshoot layer for a long 
time and some even stay there for the complete duration 
of the calculation (the dash-dotted line in Figure |3J) and 
contribute a large number of small turnover times. These 
particles, however, have only a small effect on the estimate 
of i t o for the nonrotating and slowly rotating cases where 
the turnover time might be slightly underestimated. 

From Figure [21 and the second column of Table we 
can determine that whilst l c decreases as a function of 
rotation, a similar trend is visible in u t . Thus we would 
expect that the ratio Ut/l c not to vary much as a function 
of Co. This is indeed the trend that we observe for t to 
which seems to support the conjecture that tto oc l c /u%. 

3.3.1. Dependence on the Rayleigh number 

We have made a set of runs varying the Rayleigh num- 
ber from 1.25 • 10 5 to 10 6 , denoted by ICoO, CoO, mCoO, 
and hCoO in Table |2] With the test particle method t to 
decreases slightly as a function of Ra which can be ex- 
plained by the fact that contributions of smaller scales 
are now resolved in the model leading to a shorter corre- 
lation length, indications of which are also visible in the 
velocity autocorrelations. However, we do not find any sig- 
nificant trend in the correlation time as function of Ra. It 
follows that with the present definition the Strouhal num- 
ber saturates a value approximately 0.4. This behaviour 
is similar to the Reynolds number independence of the 
Strouhal number found by Brandenburg et al. |2004). 

3.3.2. Effects of rotation 

In Paper I, we have shown that the efficiency of convection 
and overshooting is reduced as rotation increases. The lat- 
ter effect indirectly indicates that the spatial scale of con- 
vection decreases, which could lead to shorter turnover 
time if the overall velocities are less affected. The decreas- 
ing horizontal scale of convection as a function of rotation 
can be seen for example in Figure 4 of Paper I. In or- 
der to study the effects of rotation on the Strouhal num- 
ber we have made three sets of calculations with varying 
rotational influence. We denote these runs by the prefix 
CoOl, Col, and ColO, which correspond to the approx- 
imate Coriolis numbers 0.1, 1, and 10 in the runs (the 
actual values of Co achieved in the calculations vary, see 
Table We also probe the latitudinal dependence by 
making calculations at latitudes 0°, —30°, —60°, and —90° 



for each rotation rate. The correlation and turnover times, 
and the resulting Strouhal numbers for all the calculations 
are presented in Table 

For the slowest rotation, we find little differences to the 
nonrotating case CoO which was discussed above. Also the 
variation of the timescales as function of latitude is small 
in comparison to the variation within individual runs. For 
the Col set the correlation time is markedly shorter (about 
6-7 time units as opposed to 9-10 in the run CoO and 
the CoOl set). This can be explained by the deflection of 
the vertical flows by the Coriolis force, leading to smaller 
vertical scale of convection which tends to shorten the 
correlation time. Support for this conjecture is given by 
the decrease of the turnover time from the test particles. 
A similar decrease in correlation time is noted to occur 
also if it is estimated from the horizontal velocities. 

For the most rapidly rotating case, ColO, the trend 
of decreasing spatial scales continues. This is manifested 
in the clearly shorter turnover time i to as opposed to the 
more slowly rotating runs discussed above. Here, it is also 
important to note the misleading value of the turnover 
time given by the simple method, Eq. (0J. Whereas t[„ 
increases with rotation due to the smaller velocities in 
general, the actual turnover time decreases due to the 
smaller spatial scale of convection. As for the correlation 
time, the strong Coriolis forces tend to disrupt the cellular 
structure of the convection rapidly resulting in a shorter 
r c . Considering the Strouhal number, the decrease in r c 
overweights the decrease of the turnover time so that the 
actual value of St decreases as well. 

We find that St is more or less consistent with a pow- 
erlaw St oc Co -0 ' 32 as a function of the Coriolis number 
for moderate and rapid rotation (see Figure EJ). 

3.4. Implications for first order smoothing 

The relatively large values of St raise question of the 
validity of the first order smoothing approximation. To 
our knowledge, this question has been considered only by 
NS88, who derive a fourth order correlation approxima- 
tion from which they derive the transport coefficient re- 
sponsible for the generation (a-effect) and diffusion of the 
mean magnetic field for homogenuous and isotropic tur- 
bulence for which the velocity field and the correlation 
time is known. Their results indicate that the fourth or- 
der contributions, become comparable to the second order 
effects, i.e. the FOSA result, if the Strouhal number is of 
the order 0.5 or larger and that the cumulative expansion 
fails to converge if St > 1. In this light our results would 
indicate that one should be cautious in applying FOSA in 
the interpretation of the present convection calculations. 
Direct testing of this, however, requires that we compare 
analytical mean-field expressions of the relevant transport 
coefficients to the numerical results. Preliminary results 
indicate that the correlation time needed to fit the nu- 
merical data is significantly shorter than that measured 
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0.1 1.0 10 0.1 1.0 10 

Co Co 



Fig. 5. The correlation time r c (left) and the Strouhal number St (right) functions of the Coriolis number. Powerlaws 
r c oc Co - and St oc Co - are plotted. The stars, diamonds, triangles, and squares represents the calculations at 
latitudes 0° (equator), —30°, —60°, and —90° (south pole), respectively. 



from the velocity autocorrelation function (see Kapyla et 
al. l2UU5|) . 

4. Conclusions 

We estimate the nondimensional measure of the correla- 
tion time, the Strouhal number, from numerical models 
of convection. We calculate the correlation and turnover 
times separately from the velocity autocorrelation func- 
tion and the trajectories of embedded test particles, re- 
spectively, and define the Strouhal number as the ratio of 
the two. 

The Strouhal number arises in the mean-field theo- 
ries of angular momentum transport and hydromagnetic 
dynamos where its value determines the validity of cer- 
tain widely used approximations, such as the first order 
smoothing. These approximations are based on a cumu- 
lative scries expansion of the relevant turbulent correla- 
tion, e.g. the electromotive force in the dynamo theory. 
Essentially, the higher order terms in this expansion are 
proportional to the Strouhal number. Thus the value of 
St determines whether or not the expansion converges. 

The main results can be summarised as follows: 

— We find that the correlation time does not depend on 
the Rayleigh number for the parameter range explored 
in the present study. Note that the results for r c pre- 
sented here differ significantly from those of Kapyla et 
al. 1)2005(1 due to an error in the analysis in the latter 
study. 

— As function of rotation, r c decreases by a factor of 
roughly three when the Coriolis number is increased 
from to 10. At the same time, the turnover time de- 
creases only by about 20%, leading to an approximate 
relation St oc Co -0 ' 32 for moderate and rapid rotation. 

— A noteworthy fact is that if one takes the usual esti- 
mate of dividing the characteristic length by the char- 



acteristic velocity, the turnover time increases as func- 
tion of rotation due to the overall decreasing velocities. 
However, the test particle data gives an opposite result. 
The discrepancy is due to the fact that even though 
velocities decrease in general, the spatial scale of con- 
vection decreases even more. Thus it is imperative that 
the definition of St takes into account the spatial scale 
of the motions, a fact which was also noted by NS88. 

Our inability to find any dependence on the Rayleigh 
number and the relatively high values (0.1 . . .0.4) of the 
Strouhal number raise the question of the validity of 
FOSA. In the study of NS88, although not from a directly 
comparable context, the higher order effects became vis- 
ible when St was of similar magnitude. At the moment, 
we cannot draw a firm conclusion whether or not it is 
safe to use FOSA or not. We think that the next step 
in the direction of determining this would be to develop 
mean-field expressions for the transport coefficients that 
are directly applicable to the present numerical calcula- 
tions. Furthermore, these expressions should be derived 
using FOSA and preferably also with some higher order 
expansion in order to determine which of them gives a 
better description of the numerical result, and, more im- 
portantly, whether the whole concept of the cumulative 
expansion is applicable. This study, however, does not fit 
in the scope of the present paper. 
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